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Abstract 

Qh' A family of spin-lattice models are derived as convergent finite 

Q ■ dimensional approximations to the rest frame kinetic energy of a 

barotropic fluid coupled to a massive rotating sphere. In not fixing the 
angular momentum of the fluid component, there is no Hamiltonian 
equations of motion of the fluid component of the coupled system. 
• th . This family is used to formulate a statistical equilibrium model for 

rS | the energy - relative enstrophy theory of the coupled barotropic fluid 

- rotating sphere system, known as the spherical model, which because 
of its microcanonical constraint on relative enstrophy, does not have 
the low temperature defect of the classical energy - enstrophy theory. 
This approach differs from previous works and through the quantum - 
classical mapping between quantum field theory in spatial dimension 
d and classical statistical mechanics in dimension d+1, provides a new 
example of Feynman's generalization of the Least Action Principle to 
problems that do not have a standard Lagrangian or Hamiltonian. 
A simple mean field theory for this statistical equlibrium model is 
formulated and solved, providing precise conditions on the planetary 
spin and relative enstrophy in order for phase transitions to occur 
at positive and negative critical temperatures, T + and T_. When the 
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planetary spin is relatively small, there is a single phase transition at 
T_ < 0, from a preferred mixed vorticity state v = m for all positive 
temperatures and T < T_ to an ordered pro-rotating (west to east) 
flow state v = n u for T_ < T < 0. When the planetary spin is rela- 
tively large, there is an additional phase transition at T + > from a 
preferred mixed state v = m above T + to an ordered counter-rotating 
flow state v = rid for T < T + . A detailed comparison is made between 
the results of the mean field theory and the results of Monte-Carlo 
simulations, dynamic numerical simulations and variational theory. 
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1 Introduction 



Consider the system consisting of a rotating high density rigid sphere of 
radius R, enveloped by a thin shell of non-divergent barotropic fluid. A com- 
parison with the divergent case - coupled shallow water model - will be given 
in the appendix. The barotropic flow is assumed to be inviscid, apart from an 
ability to exchange angular momentum and energy with the infinitely mas- 
sive solid sphere. In addition we assume that the fluid is in radiation balance 
and there is no net energy gain or loss from insolation. This provides a 
crude model of the complex planet - atmosphere interactions, including the 
enigmatic torque mechanism responsible for the phenomenon of atmospheric 
super-rotation - one of the main applications motivating this work. We con- 
struct a simple mean field theory for the equilibrium statistical mechanics 
of this problem with the specific purpose of investigating its critical phe- 
nomenology - phase transitions that are dependent on a few key parameters 
in the problem. The adjective simple refers to the fact that in this mean 
field theory, a comparison of the free energy between two special macrostates 
in the problem - a disordered vorticity state and a unmixed vorticity state 
representing rigidly rotating flows - is made to determine the critical tem- 
peratures. Comparisonsof the results in this paper with those obtained by 
a zero temperature variational method and also those predicted by a more 
sophisticated mean field theory based on the Bragg method underscore the 
fact that this simple mean field theory is nonetheless quite powerful in pre- 
dicting the critical pehenomenology of the coupled barotropic flow - rotating 
solid sphee system. 

This paper will therefore contain new results - the simple mean field the- 
ory and its implications - together with a fairly extensive survey of the main 
ideas, approaches and results of applying statistical mechanics to complex 
geophysical flows in the past two decades. The reader is invited to consult 
the book j^T] for the details of many topics mentioned here. 

For a geophysical flow problem concerning super-rotation on a spherical 
surface there is little doubt that one of the key parameters is angular mo- 
mentum of the fluid. In principle, the total angular momentum of the fluid 
and solid sphere is a conserved quantity but by taking the sphere to have 
infinite mass, the active part of the model is just the fluid which relaxes by 
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exchanging angular momentum with an infinite reservoir. 

It is also clear that a quasi - 2d geophysical relaxation problem such as 
this one will involve energy and enstrophy. The total or rest frame energy 
of the fluid and sphere is conserved. Since we have assumed the mass of 
the solid sphere to be infinite, we need only keep track of the kinetic energy 
of the barotropic fluid - in the non-divergent case, there is no gravitational 
potential energy in the fluid since it has uniform thickness and density, and 
its upper surface is a rigid lid. 

At this point there are in principle two distinct infinite reservoirs in the 
relaxation problem, namely, the energy and angular momentum ones. The 
author's approach to this class of geophysical problems differs from previous 
works in two major aspects, namely (A) angular momentum of the fluid is 
not conserved but rather the fluid relaxes by exchanging energy and angular 
momentum with the solid sphere, and (B) the energy and angular momentum 
reservoirs is combined into a single one. (A) is justified by one of the aims of 
this approach - to model and give a relaxational explanation of the enigmatic 
super-rotation problem in the Venusian atmosphere - and by the physics 
of angular momentum of the coupled barotropic fluid - sphere systems on 
which the approach is based. Indeed, previous works on barotropic flows 
on a rotating sphere with trivial topography [3], [T3| that are based on the 
classical energy - enstrophy theory for the barotropic vorticity equation - 
hence, conserving the fluid's angular momentum - have failed to discover any 
phase transitions within the valid temperature range of the corresponding 
Gaussian models. 

The justification of (B) takes the form of the fact - easily seen by a 
calculation on the reader's part - that the rest frame kinetic energy of the 
non-divergent barotropic fluid has two distinct parts, the second of which 
is proportional to the fluid's changing net angular momentum relative to a 
frame that is rotating at the fixed angular velocity of the sphere. 

Any statistical mechanics theory on a subject matter in other than the 
traditional domain of microscopic elements, requires some discussion of what 
we mean by temperature and entropy. Statistical temperature for macro- 
scopic flows is by now a recognized scientiifc concept. Unlike the standard 
notion of temperature which measures the average kinetic of molecular mo- 
tions, macroscopic flow temperature is a measure of the average kinetic en- 
ergy contained in eddies which can vary over a wide range of length scales. 
Given a macroscopic flow state or macrostate vorticity distribution, there 
belongs a suitable flow temperature that depends on the average energy of 
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the eddies in the flow. 

Conservation of relative enstrophy is treated here as a micro canonical 
constraint, modifying the classical energy-enstrophy theories [TT], [3], [T^j . 
j^Zj in substantial ways, chief amongst them being removal of the Gaus- 
sian low temperature defect while retaining the possibility for exact solution 
of the model. Since the classical energy-enstrophy theories are all doubly 
canonical in the energy and enstrophy, they are essentially equivalent to 
Gaussian models which is the simplest of the few exactly-solvable statistical 
mechanics models known to man. Hence, previous applications of equilib- 
rium statistical mechanics to geophysical flows j^j, [Hj, |37j have largely 
used it. However, replacing the canonical constraint on enstrophy by a mi- 
crocanonical one yields a significant benefit, namely, statistical equilibrium 
models, known collectively as the spherical model, that are well-defined for 
all positive and negative temperatures. For the aims of this paper, which is 
to investigate the precise conditions for phase transitions in the equilibrium 
statistics of the coupled barotropic fluid - rotating solid sphere system, this is 
important because a phase transition could in principle, occur at any positive 
or negative temperature. 

Higher vorticity moments are considered to be less significant than en- 
strophy in statistical equilibrium models of quasi-2d geophysical flows. We 
discuss this point in greater detail in 

2 Summary of objectives, approach and re- 
sults 

It is difficult to solve the spherical model for the coupled barotropic fluid - 
rotating solid sphere system analytically in closed form. Progress on this cur- 
rent topic is discussed in and on the author's webpage: www.rpi.edu/~limc 
Other more tenable methods - both numerical and analytical ones - will there- 
fore have to be used to address this problem. This research programme will 
use three related methods to investigate the statistical equilibrium proper- 
ties of the energy-enstrophy theory of the coupled barotropic flows system, 
namely (i) a simple mean field approach - presented here - which is based on 
the notion of conserving an averaged relative enstrophy, (ii) Monte - Carlo 
simulations of the spherical model (with the microcanonical constraint on 
relative enstrophy), and (iii) Bragg mean field method. In this paper we 
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present the first. The second is presented in another paper ^H] with prelim- 
inary results announced in the AIAA paper [40J. The third is based on an 
intermediate mean field method JU] which gives a renormalized expression 
for the free energy in terms of the coarse-grained barotropic non-divergent 
vorticity, without using a relative enstrophy constraint These approx- 
imate expressions for the free energy yield values for positive and negative 
critical temperatures for the coupled barotropic fluid - rotating solid sphere 
system which are consistent with the mean field results in this paper, and the 
Monte-carlo simulations in JE] . Similar to the Curie- Weiss theory for phase 
transitions in ferromagnets ^U], the special form of the hyperbolic tangent 
plays a key role in [T7j and will be taken up again in a paper on the associated 
nonlinear fixed point result. 

The specific aim for formulating a simple mean field theory in this paper 
is to find precise conditions on planetary spin, enstrophy and energy for phase 
transitions in the equilibrium statistical mechanics of the coupled barotropic 
fluid - rotating solid sphere system. The formulation of this simple mean field 
theory is based on common features of several recent mean field theories for 
the 2D Euler equations but not equivalent to any single one of them - see 
Joyce and Montgomery [33], Lundgren and Pointin [33], Miller |44j . Robert 
051, Caglioti [HE]. 

Moreover, the above mean field theories are all known to be asymptoti- 
cally exact in a nonextensive continuum limit; proofs of this property were 
mainly constructed via large deviations methods and result in nonlinear el- 
liptic mean field equations as in [33], j3S]. Although a rigorous proof for this 
property has not yet been constructed in the case of the coupled barotropic 
fluid - rotating solid sphere system, we are certain that only the details differ 
from the general line of approach used in [33], j3H]. Our point of departure, 
then, is to assume that the mean field is asymptotically exact for the coupled 
barotropic fluid - rotating solid sphere system. Asymptotic exactness [3%] . 
[33], |33j . [33] of the mean field is a useful property for it implies that (a) a 
nonextensive continuum/ thermodynamic limit is well defined for the fam- 
ily of finite size spin-lattice Hamiltonians H^, that we will derive, as finite 
dimensional approximations of the rest frame barotropic kinetic energy of 
nondivergent flows, and (b) the thermal properties of this family of statisti- 
cal equilibrium spin-lattice models, in the thermodynamic limit as mesh size 
N tends to oo, is completely determined by the mean field. 

The mean field results in this paper are obtained without solving any 
nonlinear elliptic PDEs, an important point which will be discussed further, 
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in view of our main aim of calculating precise critical temperatures. In 
particular, using this simple mean field theory, we predict at least one and as 
many as two critical temperatures for the equilibrium statistical mechanics 
of barotropic flows on a sphere, depending on the planetary spin Q and 
relative enstrophy Q r . For relatively low values of the planetary spin Q > 0, 
there is no positive temperature phase transition, just a negative critical 
temperature T~(fl, Q r ) < between a mixed vorticity state v = m (without 
any long range order) for T, positive or negative, less hot than T~ and 
the prograde solid-body flow state v = n u for negative T hotter than T~. 
For large enough planetary spins, there is in addition a phase transition at 
positive critical temperature T c > between a hot mixed state v = m and 
the cold retrograde solid-body flow state v = n d . 

The history of applying statistical equilibrium methods to 2D turbulence 
is characterized by several distinctive findings, such as Onsager's discovery of 
negative temperatures when the vortices are confined, the presence of coher- 
ent structures and an inverse cascade of energy to large scales [TTJ , [T3] . Two 
dimensional vortex statistics also have the dubious distinction of not sup- 
porting any phase transitions at positive temperatures (cf. [33] i [12]; [TT] . 
[55j . [H] . [13] . [52] . [23]). Until recently, quasi 2D turbulence in the context 
of barotropic rotating flows over trivial topography appeared to retain many 
of these properties, including the lack of phase transition. Nontrivial topog- 
raphy appears to have qualitatively different and anisotropic effects. In their 
classic paper [37], Salmon et al derived the statistical equilibrium states of 
an inviscid unforced single-layer of fluid in a periodic box with nontrivial 
bottom topography in the beta plane approximation, and reported correla- 
tions between bottom topography and expected stream function when flow 
energy is low. Frederiksen [I-D reported statistical equilibrium results for the 
coupled barotropic fluid - rotating solid sphere system with nontrivial topog- 
raphy, that are consistent with several numerical and experimental findings, 
namely, topographic effects are asymmetric in the sense that eastward solid- 
body flows over nontrivial topography is less stable compared with westward 
solid-body flows. 

However, when the topography is trivial, Frederiksen and Sawford [3] 
reported that the eastward and westward solid-body flows are both stable 
but did not find any phase transitions. By going to a simple but powerful 
mean field approach, not based on solving PDEs numerically, we calculate 
asymmetric positive and negative critical temperatures in terms of planetary 
spin and relative enstrophy. We believe that our result on the existence of 
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a positive critical temperature in the rotating barotropic flows system when 
the planetary spin is sufficiently large, and a negative critical temperature 
for all values of the planetary spin, is the first in the large literature on 
the equilibrium statistical mechanics of geophysical flows since the papers 
by Salmon, Holloway and Hendershott [37], Frederiksen and Sawford [3] and 
Carnevale and Frederiksen |39|. 

This paper is organized into the following main sections and appendix: 
(1) the extended Planck's theorem, (3) derivation of Ising model-like spin- 
lattice Hamiltonian for the kinetic energy of the barotropic flows, (4) the 
classical energy-enstrophy theory based on this spin-lattice model, (5) the 
spherical model for the barotropic flows on a rotating sphere and why it is 
not exactly-solvable, (6) simple mean field theories of the spin-lattice Hamil- 
tonian and its thermal properties, (7) differences between the Frederiksen 
and Sawford theory [3] and the spherical model, (8) conclusion and (9) ap- 
pendix: background and properties of the coupled barotropic fluid - rotating 
solid sphere system, 

3 Equilibrium statistical mechanics 

Planck's theorem will be used to find the most probable equilibrium states 
in an isothermal setting: minima of a Gibbs free energy per lattice site 
/ determines the thermodynamically stable states. It is wellknown since 
Onsager's seminal paper |35] that 2D vortex statistics is characterized by 
negative temperatures at high values of flow kinetic energy. Since we are 
interested in exploring the entire range of flow kinetic energy in our model, 
we will have to extend Planck's theorem for negative temperatures T: 

Extended Planck's Theorem: The most probable (thermodynamically 
stable) state at negative (positive) temperature T corresponds to maxima 
(minima) of the Gibbs free energy per site 



f = U-T3 



where 



U=(H) 



is the internal energy per site and 



S 
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is the mixing entropy per site, given in terms of the probability distribution 
p(s) for site value s. Here ks is the Boltzmann constant. 

We sketch a proof next, using an alternative definition of the total entropy 
due to Boltzmann. The total free energy F = U — TS is related to the Gibbs 
partition function by 

Z(P) = e ~' 9ir(m) = e -P(U{m)-TS{m)) ^ 

where 

F(m) = U(m) — TS(m) 
denotes extremal values of the free energy, 



is the inverse temperature and 

S(m) = fcg In. W(m) 

is the total entropy in terms of the degeneracy (or number of microstates) 
W(m) of the most probable state v = m. Then 

Z(j3) = e -W(m)+lnW(m) = W ( m ) e ~/^M = p( m ) ( 2 ) 

is equal to the Gibbs probability of the most probable state m. From (0) 
and (J2J), we deduce that maximal (minimal) values F(m) of the total free 
energy give the dominant contribution to Z(f3) for negative (resp. positive) 
temperature. 

Remark 0: From the proof of the above theorem, it follows that the 
most probable state, also the thermodynamically stable state, is not necessarily 
the one, call it m' that maximizes W = exp (S/ks) , but rather that which 
maximizes the product W / exp(— f3U). Nonetheless, at temperatures T where 
\T\ is small, it is likely that the state extremizing F is closer to one, call it m 
that extremizes U instead of m! maximizing S; and when \T\ is large, it is 
likely that the state extremizing F is closer to m' that maximizes S than that 
extremizing U. We note the significant point that for systems that support 
negative temperatures such as barotropic flows in particular, the specific heat 
is nonetheless positive even when the temperature is negative. 
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A precise correspondence exists between the statistical equilibrium prop- 
erties of barotropic flows on one hand and its dynamical properties on the 
other. This correspondence is encapsulated in the Minimum Enstrophy or 
its equivalent, the extremal kinetic energy principle, which form a pair of 
dual variational principles for the steady-states of the barotropic flows (cf. 
Leith j2H] 5 Young [2j, Prieto and Schubert j2H] and Lim [151). The maxi- 
mal kinetic energy steady-state of the barotropic flows is closely related to 
the most probable flow state at hot enough negative temperatures in the 
corresponding statistical mechanics model. Likewise, the minimal kinetic 
energy steady-state of the barotropic flows is related to the most probable 
flow state at low positive temperatures. At temperatures sufficiently close to 
0, this correspondence follows directly from the form of the Gibbs free energy 
/ = U — TS ~ U, no matter which of the two general categories below, the 
barotropic flows thermal system happens to be in. 

In general, the Gibbs canonical ensemble consists of the standard form 
for the probability measure 

PgW = -^^exp(-/3#M -iiT[w}) 

where H[w] is the energy and (3 an inverse temperature while T[w] represents 
a significant conserved quantity and \i a chemical potential or Lagrange Mul- 
tiplier conjugate to T[w]. The partition function or configurational integral 
Z[f3, fj] provides the normalization required to make -PgM a probability mea- 
sure. It is customary to include only the key conserved quantities as canonical 
constraints T[w]. Indeed, in the case of quasi-2D turbulence, there are an in- 
finite number of conserved quantities, namely the moments of the relative 
vorticity or higher order enstrophies, and it is unnecessary to include all of 
them as canonical constraints in the Gibbs probability Pq. 



4 Spin-Lattice Approximation 

Given the well known fact that Gibbs' canonical ensemble and the corre- 
sponding partition function for the spherical model - to be discussed in detail 
below - are closely related to path-integrals and therefore extremely complex 
mathematical objects, a rational approximation procedure based on finite 
dimensional spin-lattice models or something similar, will have to be devised 
to simulate their critical phenomenology on the computer as well as to solve 
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them by analytical means. Such a rational approximation scheme must sat- 
isfy two basic requirements when the size or order of the approximation is 
taken to infinity: (A) the resulting family of finite dimensional models con- 
verge to the correct energy functional and constraints of the problem and (B) 
the thermodynamic limit - in this case, the nonextensive continuum limit - 
of this family of approximate models exists. (A) is shown to be true for the 
family of spin-lattice models given next. (B) turns out to be true if exact 
solutions to the spherical models - to be obtained by the Kac-Berlin method 
of steepest descent in £17] - yield valid free energy expressions in terms of 
the associated saddle points in the nonextensive continuum limit. For the 
purpose of this paper, the assumption of the validity of (B) is subsumed un- 
der the earlier assumption that the mean field is asymptotically exact in this 
class of problems for the coupled barotropic flows on a rotating sphere. 

Using a uniform mesh M of N points {x\, Xjy} on S 2 and the Voronoi 
cells based on this mesh - see Lim and Nebus [20J - we approximate the 
relative vorticity by 

JV 

w {x) = ^2sjHj(x) 
j=i 

where 

Sj = w(xj) G (—00, 00) 

and Hj is the indicator function on the Voronoi cell Dj centered at Xj, that 
is, 

Hj(x) = 1 for x e Dj, 
Hj{x) = for x i D r 

The real valued spins Sj should henceforth be viewed as coarse-grained 
or block averaged vorticity - resulting from a single step renormalization 
procedure outlined in [T7j - that is ideally suited to the mean field approach 
in this paper. Another point of view describing the real valued spin states 
{sj} as macrostates can be found in the book jH]. 

Recall that the Voronoi cell Dj is defined to be all the points in S 2 that 
are nearer to Xj than to any other points in the mesh M. This confers the 
essential property that the Voronoi cells is a disjoint cover for S 2 , that is, 

DjHD k = 0, 
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N 

u 

J'=l 



Uniformity of the mesh M confers the other essential property that the areas 
Aj of the cells are equal, that is, 

A 3 = \D 3 \ = ^. 

N 

The rest frame kinetic energy given in terms of the parameters of a frame 
that is rotating at the fixed angular velocity of the solid sphere has been 
shown to be 



H 



1 



dx (w + 2Q cos 9) G(w) 

where the fundamental solution of the Laplace-Beltrami operator on S 2 is 

tp(x) = G(w) — dx In |1 — x • x'\ w(x'). 

Under the above approximation for the relative vorticity w, the truncated or 
lattice approximate energy 



H N = 



\ J s2 dx (^2s J H J (x) + 2ncos9 ] j G (jtsMx)^ 



N N r 

EE 



S 2 



2 ^ ^ 

Z j=l k=l 

H as N — > oo 



dx Hj(x)G(H k ) 



N r r i 

SjS k — fij^ / dx cos 9G (Hj) (3) 
, =1 US' J 



under the calculus rules of (Lebesque) integration. With the interactions 

J jk = [ dx Hj(x)G (H k ) 



and the external fields 



Fj = Q f dx cos9G(H j ), 

J 
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the truncated energy takes the standard form of a spin lattice model Hamil- 
tonian, 

I N N N 

ff W = -nEE JjkSjSk - F 3 S 3- ( 4 ) 
3=1 k=l 3=1 

The interactions Jjk are logarithmic and thus, long range. The external fields 
Fj are non-uniform and linear in the spin Q > 0, and represent the coupling 
between the local relative vorticity or spin Sj and the planetary vorticity field 
pcos9. They are turned off for the special case of single layer inviscid vortex 
dynamics on a non-rotating sphere. The presence of this inhomogeneous 
term when Q > 0, is the source of the much richer mathematical and physical 
properties of the coupled barotropic flows on a rotating sphere. 

Evaluating the integral and using the essential properties of Hj we obtain 

Jik = / dx HAx) I dx' In 11 — x ■ x'\H k (x') 
Js 2 Js 2 

167T 2 , | . , 

— > m |1 — Xj ■ Xk\ as iv — > oo. (5) 

For the integral in the external fields Fj we obtain 

Fj = £1 J ^dx cos# J 2 dx' In |1 — x • x'\Hj{x') 

= / dx'HAx' ) I dx cos ^ In 11 — x ■ x'\ 

Js 2 Js 2 

after using the symmetry of the inverse G to the Laplace-Beltrami operator 
on S 2 . Then 

Fj = Oncost / dx' H] {x')j d.rrj.n\, , 

= ~-Q\\cos9\\ 2 I dx'Hj(x')ilJ l0 (x') 
2 J s 2 

2n 

— > — — Q\ | cos 6\ 1 2^10 (^j) as N — > oo, (6) 
where || cos 0| | 2 is the L 2 norm of the function COS6 1 , that is, 



cos#|| 2 = J J dx cos 2 9, 
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and the spherical harmonic ip 10 which represents the relative vorticity of 
solid-body rotation, satisfies the eigenvalue relation, 

|| cos 6* 1 1 2 

\ lm = -1(1 + 1) = _2 for I = l,m = 0. 
The truncated relative enstrophy is given by 



N 



dx w 2 
s 2 



J s2 dx \^ s i H A x )\ 

-^5>J-r r asiV-oo, (8) 



4tt N 



after using the two essential properties of the indicator functions Hj. 
Lastly, the truncated total circulation is given by 

JV 



TCn = J ^dx w = J 2 dx SjHj(x) 
4tt n 

= — Sj — > TC as N — ► oo. 

5 Classical and Recent Energy-Enstrophy The- 
ories 

To put the mean field theory here in perspective, we briefly review the two 
main theories - both of which are based on Kac's inventions circa 1952 - used 
to investigate the statistical relationships between energy and enstrophy in 
macroscopic flows. The first is the Gaussian model which forms the basis 
of all the classical energy-enstrophy theories. The second is the spherical 
model for which the canonical enstrophy constraint is replaced by the micro- 
canonical form. Versions of the spherical models for new energy - enstrophy 
-circulation models of macroscopic flows were first introduced by the author 
beginning in early 2000 - see the new Springer- Verlag book by Lim and Nebus 
which will be publsihed in October 2006 |41j . 
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5.1 Gaussian model 



Almost all the early papers on variations of this model used a spectral formu- 
lation in terms of a truncated set of orthnormal eigenf unctions for the flow 
domain. Many authors presented results on the applications of this theory 
to a large range of topics in geophysical flows including two-layer flows over 
nontrivial bottom topography, and quasi-gestrophic f-plane and beta plane 
flows [37] . We will use a spatial lattice formulation instead for reasons already 
mentioned in the introduction. It is easy to show that this classical energy- 
enstropy theory is identical to the well-known Gaussian Model introduced by 
Kac which is exactly-solvable. 

The classical energy-enstrophy theory as formulated using a spatial dis- 
cretization is now given in terms of the truncated energy fl2J) and relative 
enstrophy (JHJ) by the Gaussian partition function, 



JN ~ w 



\ 1 r I N \ 

J [n ds 3 ] exp -b £ s) exp (-/3H N [S; fi]) , 



b\ N/2 r ( v 

^ I [n ex p ( - i 

in terms of the spin (vorticity) state 

S = {si,....,s N }, 
Sj G (—00,00). 

In this model, the standard deviation of the above Gaussian is given in terms 
of the parameter 

N 

where /1 is the chemical potential or Lagrange multiplier associated with 
the relative enstrophy constraint. The factor (^) is needed to make 

W[S] = ^ exp (— bY^j = i sf) a probability distribution. 

The exact solution of this Gaussian model where it is well-defined, anchors 
many of the previous works on statistical equilibrium in geophysical flows 
(cf. Salmon, Holloway and Hendershott [37j . Frederiksen and Sawford [3] 
and Carnevale and Frederiksen [39J, amongst many others). 

It turns out as will be shown next, that the choice of a canonical constraint 
on total kinetic energy H[q] and a microcanonical constraint on the relative 
enstrophy T r gives a spherical sodel formulation which although difficult 
to solve analytically, is amenable to numerical simulations and mean field 
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methods. More important is the fact that the spherical model is well-defined 
for all temperatures positive or negative, while the Gaussian model is not 
defined for certain temperatures. 

5.2 Spherical model for the coupled barotropic flows 

We need to fix the low temperature defect of the classical energy-enstrophy 
theory of the barotropic flows which was mentioned above. We will do this by 
replacing the canonical constraint on relative enstrophy by a microcanonical 
constraint, which yields a version of Kac's spherical model for the spin-lattice 
Hamiltonian (JUJ). The spherical model formulation of an equilibrium 
statistical energy-enstrophy theory for the barotropic flows on a rotating 
sphere is based on the spin lattice partition function, 

Z N = J (Ud Sj )6(NQ r - J2 4) ^p [—f3H N ] (9) 

where H^, Jjk and the external fields Fj are given in 0, © and ©• 

The microcanonical relative enstrophy constraint takes the Laplace inte- 
gral form and gives 

I r ra+ioc 

Zn = / (Hdsj) / dr/exp 

Ani J Ja-ioo 

in terms of the matrix 

K jk = r]- (3 J jk 

and the spin vector S = {s 1 , s^}. 

6 Mean Field Theory 

The spherical model for the energy-enstrophy theory of the rotating barotropic 
flows is difficult to solve but it fixes the above foundational difficulty of the 
classical energy-enstrophy theory of the coupled barotropic flows. We there- 
fore use mean field methods to derive thermal properties of the spin-lattice 
Hamiltonians, 

j JV JV N 

H N = ~« J jkSjSk - F 3 S J- 

Z j=l k=l j=l 



±qNQ r -±(S\K\S) 





N 


exp 
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This yields physically significant phase transitions for the new energy-relative 
enstrophy theory of the coupled barotropic flows on a rotating sphere. The 
main step of this mean field theory is to calculate the change in free energy per 
site between two fundamentally different vorticity states, namely, the mixed 
macrostate where the coarse-grained vorticity of opposite sign are juxtaposed 
in such a way that a site has equal probability of having spin ±s and the 
unmixed macrostate where the coarse-grained vorticity of opposite sign are 
separated into hemispheres on S 2 in a way that a site has probability (resp. 
probability 1) of having spin s (resp. — so) - depending on where the given 
site is located relative to the hemispherical pattern. 

Many detailed aspects of the formulation below can be changed without 
any loss of accuracy or precision to the stated results. Some of these are the 
neighborhood size z which determines the range of the interactions in the 
model, and the simplified mean spin values ±so- For the energy-enstrophy 
spin-lattice models, the actual neighborhood N(j) equals the set of all other 
sites in the mesh. However, although this signify a global interaction in H N , 
it is a short range one by the accepted definition - the integral of interactions 
over all sites is finite - and therefore, its mean field properties are not depen- 
dent on z. Choosing the simplified mean spin values to be ±so appears to 
be a more serious restriction - changing the real valued spins in the original 
lattice models H N to binary spins that make the resulting model more of an 
Ising type. Again, it is not hard to see - by a single step of spatial renormal- 
ization or block spin averaging - that the binary spin models we use in the 
ensuing mean field theory are equivalent to the original. 

6.1 Set-up for coupled barotropic flows 

Since the coarse-grained spin or vorticity s (x) satisfies the following integral 
relations, 



we assume, for the mean field theory, the existence of a site dependent volume 
fraction or probability distribution function v x (sq) such that 





M 



(10) 
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where M denotes the limiting values of the coarse-grained spin Sq(x) over S 2 
and 



p pm 

I dx j ds u x (s )sl = Q (11) 

JS 2 J-M 

p pM 

dx ds v x (s )s = 0. (12) 

Js 2 J-M 

Next we assume that the mixed state and unmixed mean field state, 
denoted respectively by 

v — m,n, (13) 

are characterized by m where the mean spins sq{x) are independent at 
neighboring sites and n where the neighboring spins are correlated, that 
is, s (x) = s (x') G {sq }. Equivalently, m is characterized by 

z/™(±s ) = - for all x G S 2 ; (14) 
2 

and n by 

v n x+ (s Q ) = 1 for all x + G S 2 , 
u^{so) = lforallx-GS 2 

where x are lattice sites in respectively the positive and negative hemi- 
spheres of the unmixed state n. 

The notation N(j) denotes the neighborhood of site j, that is, all lattice 
sites k which are connected to site j. Let |iV(j)| = z be the common size of 
neighborhoods N(J) or coordination number for the lattice. By varying z we 
can model interactions of different range in a variety of spin-lattice models. 
Let e denote the interaction energy scale (obtained by averaging Jjk over 
N(j)) for Hjy. For the purpose of modelling the energy-enstrophy theories, 
e < which represents an anti-ferromagnetic interaction. 

In addition to the spin-lattice Hamiltonians H^, we will need a crude 
lattice approximation for the per site mixing entropy. Following the same 
approach for Jjk and Fj in (pj) and © respectively, we derive the lattice 
approximation for the total mixing entropy in the ensuing steps, 



p pM 

S[u] = —ks dx dsv x ( s) In i/Js) 

JS 2 J-M 
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N r rM 

kB^2 J dx Hj(x) J dsv x (s)\nis x (s 



E 



dsUxAs) In u x As] 



6.2 Proofs and results when the solid sphere does not 
rotate 



First we compute the entropy per site for the mixed state v — m, 



Sn 



Ank B r M . . 47rfc? 



dsv x In v x 



In 2. 



N J-M N 
Similarly, the per site entropy in the unmixed state v — n is given by 



S n = 0, 



(15) 



(16) 



because by the defining property of the unmixed state v = n, neighboring 
mean values sq(x) are correlated perfectly. 

Next we compute the mixed state internal energy per site, 



u. 





cAf 




ds 




-M 


ez 


' rM 


T 


J-M 



l 2 







(17) 



by virtue of property (j!4jl . We compute the unmixed state per site internal 
energy by using the property that in the unmixed state, neighboring mean 
spins s(xk) satisfy 

s(xk) = s(xj) = ±s for all k £ N(j) 
with \N(J)\ = z for all j = 1, N. 

Then in this mean field theory, 



ez 



2N 



r t M ez f M 

dx ds v x sf, = / dsnv x si 

JS* J-M 2 J-M ° ° 



(18) 
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after using (fTTj). 

Given the enstrophy Q > and temperature T, and treating the mean 
spin distribution v x and state v = m, n as free parameters, the isothermal 
free energy difference per site between the mixed and unmixed states is given 
by 

^■f I'm fn 

ez ^ 4nk B T , „ , . 

= —Q — In 2. (19) 

8n N V ; 



Equivalent ly, 



fm TS, 

4irk B T 



In 2; (20) 



and 



fn TS n 



= - £ £Q. (21) 

We have nearly finished proving the following mean field result for the 
non-rotating barotropic flows using the extension of Planck's theorem that for 
negative temperatures T < 0, thermodynamically stable statistical equilibria 
corresponds to maximizers of the free energy: 

Theorem 5: If e < 0, then (i) for all T > 0, the mixed state v = m is 
preferred. For e < and T < 0, there is a finite N critical temperature 

T C (N) = < (22) 

such that if (ii) 

> T > T C (N) 
then v = n is preferred, and (Hi) if 

T < T C (N) < 0, 

then v = m is preferred with maximum mean spin entropy 

S ms =^k B \n2. (23) 
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(iv) In the nonextensive continuum limit - assumed to exist- as N 
oo, T C (N) tends to a finite negative critical temperature 

T ^ = § Is 2 dxjj w (x) J s2 dx' ip 1Q (x') In |1 - x ■ x'\ 



§ Is 2 dxif> 10 (x)G (V'io) (») 



i - ^) = Q_ , 



> T > T c 

t/jen v = n is preferred and if 

T < T c < 0, 

t/ien v = m is preferred. 

Proof: The proof of (i) follows from the above calculations. To prove 
cases (ii) and (iii), we compare the per site free energy 

T 

/m(max) = - — S ms > (24) 
with the same in the unmixed state v = n, 

F Z 

/ n (max) = --Q > 0. (25) 

07T 



When 



■-S <--0 



or equivalently > T > T C (N), (ii) v = n is preferred, and vice- versa for 
(iii). Since £ < 0, Q > and by definition, S ms > 0, the critical temperature 
T C (N) < 0. 

The proof of (iv) follows from consideration of the nonextensive contin- 
uum limit under which the free energy 

ezN . Q 



8vr 



-Q — > — [ dxip 10 (x) [ dx' ip l0 (x') In |1 
2, J s 2 J s 2 



X ■ X 



Q 
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as the number N of lattice sites tend to oo, together with the fact that the 
denominator in T C (N) - total entropy S ms in - does not depend on N. 



Remark 1: When the interaction energy scale e < 0, there is a negative 
temperature transition between the mixed state at hot T < T c < and the 
unmixed state at very hot T c < T < 0. 

Remark 2: The mixed state free energy per site is entirely entropic. The 
unmixed state free energy f n is purely an internal energy term which is linear 
in the enstrophy Q. 



6.3 Mean field theory when the solid sphere rotates 

Let us denote the two parts of Hn by 



Z i=l fc=l 

( 2) 27rfi 

Then we write the corresponding two parts of the internal energies, 

u = u (l) +u (2) 



j iV JV 



with 



■it 



L 1} = 



in 



by (fTTjl and (fT%|) respectively. 

Next we assume that the mixed state v — m is defined in (II Hj) and the 
unmixed mean field states denoted by v — n, n u and are characterized by 
correlated neighboring spins, s(x) = s(x') G {s^}. Moreover, the unmixed 
states v = n, n u and are both exactly correlated by hemispheres into 
opposite values Sq such that 

4 + so = (26) 
2vr((4) 2 + (^o) 2 ) = Qr- (27) 
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Solving for Sq we get 

4 = so = ■ (28) 

The unmixed states v = n u and are further characterized by hemi- 
spherically correlated spins Sq which satisfies the expressions f)34p and (J55j) 
respectively, that is, aligned and anti-aligned with the northern and southern 
hemispheres corresponding to the planetary spin Q > 0. 

The entropy calculations are exactly the same as when planetary spin 
Q = because the entropies S m and S n depends only on the statistical 
distribution v x (so) of the mean field relative vorticity s . 

Next we compute the values of the per site internal energies u)$ and uffi 
due to the nonzero planetary spin Q > of H N : using definition (fT^j) for the 
mixed state v = m, we get 



u 'm 



2nQ /» ' 
E SjCosOj 
\j=i i 



N 2 



ft I \ 

3=1 
2ntt ( r M 



N 2 




rM \ N 

J dsu x (s) sj Y, c os9j 



since (JHJ). Denoting by 



the lattice sites and co-latitudes that fall respectively into the hemispheres 
determined by the correlated means Sq , we obtain for the generic unmixed 
state v = n, 



2ntt /A \ 2Titt N 



u, 



(2) 

iV \i=i / JV j=i 



J2 Sj cos Bj ) = E ( s j) cos 9 j 



2ttQ 
IP 

2nQ 



[E Uxp)00Bff}\+ (E (SOO)COS^ 



\3=1 

N/2 

E 





(N/2 




+ E 




\fe=l 


AT/2 




E 


COS# fc J 


fc=i 
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where the hemispheres defined by Sq need not be aligned with the northern 
and southern hemispheres corresponding to the spin Q > 0. 
Next we compute the per site free energies 



fm U m TSj 

4nk B T 



in 



N 

and 

fn T S n 



In 2 (29) 



P7 9ttO / N/2 N/2 

= ~Qr + 9? ( 4 £ C08 6i; + S Q £ 008^ ] . (30) 



Therefore, the per site change in free energy is now given by 



A/ = f m -f r 

A7ik B T 
N 



In 2 



ez 2nQ ( N/2 _ N/2 



Remark 3: Application of this formulation to the non-rotating barotropic 
flows on the sphere at both positive and negative temperatures requires con- 
sideration of the nonextensive continuum limit under which 

ez = e(N)z(N) -> C^iV" 1 ) < 0, 
2tt n/2 2fi N/2 

— £ 4cos^+ + — E So -cos^^O(7V-i). 

j=i fc=i 



6.3.1 Positive temperature 

For T > 0, we will need to compare the minimum free energy per site in the 
mixed state v — m, 

/m(min) = -^Sms= ln2 (31) 

24 



with the same in the unmixed state v = n, 

ez 2nQ ( N/2 _ N/2 \ 

/ n (min) = ~-^Qr + -jyj- 14 Y cos9 t + s oE s o cos9 k J ■ (32) 

The extreme value of / n (min) is obtained at the most negative value of 
2nQ 



N 2 



N/2 N/2 
3=1 k=l 



which occurs when the hemispheres associated with the correlated means 
are anti-correlated with those corresponding to the planetary spin Q > 0. 
For any T > 0, if 

2nQ ( N/2 _ N/2 \ ez T 

— \ S + Y COs6 t + S Y COS0 k > ^Qr ~ J^Sms (33) 
\ 3=1 k=l J 

then v — mis preferred; if the inequality is reversed then an unmixed state 
is preferred. 

Thus, if e < 0, then for all T > and Q r > 0, the RHS of JHD is negative. 
This implies that for all T > 0, Q r > 0, and Q > 0, the mixed state v = m is 
preferred over any v = n u state where the hemispheres associated with the 
means Sq are correlated with those corresponding to spin Vt > 0, i.e., 



1 



N/2 N/2 

^ [s+ Y cos % + s o E cos ej; } > 0. (34) 
ly » 3=1 k=l 



Next, using (}3*3*j) . we compare the mixed state v = m with the unmixed 
states v = rid where the hemispheres associated with the means are anti- 
correlated with those corresponding to spin Q > 0, i.e., 



1 

N 



N/2 N/2 

sj}" Y cos ®t + «oE cos °k\ < 0- (35) 

3=1 k=l 



For e < and any fixed relative enstrophy Q r > 0, there is a positive finite 
size N critical temperature depending on Q and Q r , 

T C (Q, Q r ; N) = — 8 - Vr > (36) 

•Jms r ) 
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where 



N/2 N/2 

oo < J_ = min j- \s£J2 cos9 t + s o E cos9 k ] < °> ( 37 ) 

S iV V j=l k=l 



provided the spin f2 is large enough relative to Q, r , that is, 



FZ 

Q > n + (Q r ) = ——NQ r (38) 



which converges to a finite limit as iV — > oo. The minimum in I (|37j) is 
taken over all possible orientations of the hemispheres associated with s , 
and by J2BD, 

i- = o{Jo~ r ). 

In summary, (a) for T > T C (N) > 0, the mixed state t> = m is preferred 
over any unmixed state v = n u , and (b) for < T < T C (N), the unmixed 
state v = rid satisfying (J33j) is preferred over the mixed state v = m; (b) 
in the nonextensive continuum limit as iV — ► oo, the positive finite size N 
critical temperature 

T c (0, Q r ; N) -> T c (0, Q r ) < oo 

and moreover, if 

FZ 

Q > n + (Q r ) = ^rprNQr 

holds, then 

T c (n, Q r ) > 0. 



Remark 4: Thus, there is a positive temperature phase transition if the 
planetary spin is large enough. 

6.3.2 Negative temperature 

An extension of Planck's theorem to negative temperature states that for 
T < 0, the state with maximum free energy is preferred. Thus, we compare 
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the per site free energies of the mixed v = m and unmixed states v = n u and 
n d : 



/m(max) 



T 



f n (max) = ~^Qr + ^ 



N/2 N/2 
3=1 k=l 



For T < 0, 



£2 27rfi 



T 



N/2 N/2 

4 £ cos 0+ + s - J2 cos 0* ) > -j^S ms > (39) 



implies that the unmixed states are preferred. Solving this inequality for the 
finite size N critical temperature we get 



T~(Q, Q r - N) = ±-M < o 



(40) 



where 



1 



N/2 N/2 

oo > I + = max — 4 E cos ^ + + s o E cos 1 

s " N ' j=l '' k=l 



-I- > (41) 



and 

I + = 0{^Q~ r ). 

In this case when (c) T < T c r(iV) < 0, the mixed state v — mis preferred over 
any aligned unmixed state v = n u and thus, over any anti-aligned unmixed 
state v = rid as well since the left hand side of (J3~9")) evaluated at v = rid an d 
v = n u satisfies 



ez 2-nVL 



N/2 N/2 
j=l k=l 



(n d ) < 



ez lixVl 



N/2 N/2 

cos ^ + + s oE cos ^ 

3=1 k=l 



T 



When (d) T c (N) < T < 0, the aligned unmixed state v = n u is preferred 
over v — m. It is important to note that a negative finite size N critical 
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temperature T c (Q, Q r ; N) < exists for any spin Q > 0, unlike the positive 
finite size critical temperature T C (N) which exists only for spins satisfying 
(EH). 

Remark 5: Under the nonextensive continuum limit N — > oo, 

T c 7(n,Q r ;N) -> T-(fi,Q r )<oo 
T-(Q,Q r ) < 0, 

because in MUj) , the denominator is negative by definition and does not depend 
on N, and the numerator tends to a finite positive limit by virtue of Remark 
3. This implies that there is a negative temperature phase transition for all 
values of the planetary spin. 

Since the numerator in (J4(Jj) contains the positive term — ^NQ r (for e < 
0), we can define a second negative critical temperature 

27rfi/_ - §^NQ r , 
T-(n,g P ;JV) = < (42) 

"ras ) 

provided the spin Q > is small enough, that is, 

EZ 

Q < n+(Q r ) = TTTYr-NQr. (43) 
l07T z i_ 

We note that T c (Q,Q r ;N) < and T C (N) > in (jSHJ) are the same ex- 
pression corresponding to the two sides of the equality in (JJ3J). In this case 
when (e) T < T c (JV) < 0, the mixed state v = m is preferred over the 
anti-aligned unmixed state v = n d and when (f) T c (JV) < T < 0, the 
anti-aligned unmixed state t> = is preferred over v = m. 

Comparing (j4*Uj) to (J4"2~j) . we deduce that when the spin f2 satisfies (JMJ), 

T-(JV)<T c -(JV)<0. 

By comparing the per site free energies of unmixed states v = n u and 
when T < (more specifically T~ < T < T c < and T~ < T c < T < 0), 

tut \ ez ^ , 2?rfi r 

/«W = ~^^ + ^r J + 
/»W = -— Q r + — 1_, 

we deduce at once that the aligned unmixed state v = n u is always preferred 
over the anti-aligned unmixed state v = n<i for T < 0. 
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Remark 6: Under the nonextensive continuum limit N — > oo, 



r--(Jv)->r c -<oo 



and 



T 6 7 < T c - < 
provided the planetary spin is small enough, that is, 



n < tt + (Q r ) 



ez 



NQ r . 



16n 2 L 



Collecting together the above computations, we summarize the phase 
transitions in this simple mean field theory for the coupled barotropic flows 
- rotating sphere system: 

Theorem 6: (A) for large enough spins Q > 0, the anti-aligned unmixed 
state v — rid changes into the mixed state v = m at T c (Q,Q r ) > 0, (B) 
this mixed state v = m continues to be preferred for all positive T > T c 
and negative T < T~, changing to the aligned unmixed state v = n u at 
T~{VL,Q r ) < 0, and (C) (i) for large enough spins Q, the unmixed state 
v = n u persists for all T such that T~ < T < 0, and (ii) for small enough 
spins Q < Q + (Q r ), the state v = n u persists as the preferred state for all T 
such that T~ < T < 0, but for T such that T c < T < 0, the state v = na 
is preferred over v = m but not v = n u . 

In the nonextensive continuum limit, the mean field critical temperatures 
of the energy-relative enstrophy theory for this problem: 



oo > 





min WiQr (f S 2 dx ttC^ w G[w(x')}) + 



ms 



precisely when 



> 0; 




2C 
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and for all planetary spins fi > 0, 



— oo < T c (Q, Q r ) = lim 



2nM + - §NQ r 



N-+00 (~S ms ) 

max Wj Q r (Js'zdx Qcos6(x)G[w(x')]) + % 
( S ms ^j 

-\VLC-fQ~ r - -f ^ ^ 



<7 

where J* are given by (JIT]) and (j3*7|) . C = 1/ J 52 rfx cos 2 8, the min^^^ (resp. 
max„, j Q r ) is taken over all relative vorticity w(x) with fixed relative enstrophy 
Q r and the total entropy is 

S ms = -A-Kk B \n2 > 0. 



Proof: Using the definitions of the interaction energy scale e(N) and the 
size z(N) of the interaction neighborhood, we deduce that the term —^NQ r 
in the numerator of T C (N), being the finite dimensional representation of the 
per site spin-spin interaction internal energy in the unmixed states v = n u 
or nd (as opposed to the spin-f2 interaction term), tends to 

dxi/) w (x) [ dx'ip 10 (x')\n\l-x-x'\ = ^ > 0. (44) 

Thus, from the definition of the minimum and the derivation of the 
spin-lattice Hamiltonians Hn from the rest frame pseudo kinetic energy H 
of the coupled barotropic flows, the finite N critical temperature 

un N) = 2^min^(4Efi 1 2 cos^ + + , cob**) - § NQ r 



( S ms ) 

min. ^ ( f„ n drr. O 

T c (n,Qr) 



min^Q^. (/ 5 2 dx Q cos 9(x)G[w(x')}) - ^ 



( S ms ) 

where the denominator does not depend on N. The minimum in the numer- 
ator is taken over all relative vorticity w(x) with fixed relative enstrophy Q r , 
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and is a well-defined finite negative quantity that is proportional to f2 and 

to \[Q~r - 



mm 

W,Q r 



f dx QCip w G[w(x')]\ = f dx ttCi) l0 G[J~Q r il) 

Q r dx ib 
Js 2 



^VLC\jQ r 1^ dx ip 2 w 



Since 



EZ 



_> n°°(Q r ) = — > o, 

when Q > Q°^(Q r ), the mean field critical temperature 

T c (Q,Q r ) > 0. 
A similar argument based on 

max (^J 2 dx QCtp 10 G[w(x')]j = J ^ dx nCt/j 10 G[—\jQ^^ 10 \ 
proves the existence of the continuum limit of T~ to conclude the proof. 



It is useful to restate and separate the above phase transitions into two 
categories depending on planetary spin. 



Theorem 7: When the planetary spin Q < Q + (Q r ), (i) there is no 
positive temperature phase transition and the mixed state v = m is preferred 
for all T > 0, (ii) there is a negative temperature phase transition at T~ < 
(which exists irrespective of the value of planetary spin), where the preferred 
mixed state v — m for all T < T~ < 0, changes into the aligned state v = n u 
which is preferred over both v — m and v = for all T~ < T < 0, and (Hi) 
there is a secondary transition at the hotter temperature T c < (that is, 
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T~ < T c ), where the intermediate state v = rid changes place with v = m 
in order of thermal preference. Letting -< denote 'has smaller free energy 
than', we summarize the state preference for the case Q < Q + (Q r ) : 

n-d -< n u -< m for T < T~ < 
n d -< m -< n u for T 6 7 < T < T r 
m -< rid -< n u for T c < T < 0. 

When the planetary spin Q > Q + (Q r ), (iv) there is a positive critical 
temperature T C (Q, Q r ) given by \3b\) at which the preferred state changes 
from v = rid f or 0<T<T c tov = m for all T > T c and all negative 
T < T~ < 0, and (v) there is a negative critical temperature T~ < that 
exists irrespective of the value of planetary spin, at which the preferred state 
changes from v = m to v = n u for all negative T > T~ . We summarize the 
state preference for the case Q > Q + (Q r ) : 

nd -< n u -< m for T < T 6 r < 
rid ~< m -< n u for T~ < T < 0. 

7 Conclusion 

A related approach used to derive nonlinear stability properties of steady- 
states in the coupled barotropic flow - rotating sphere system is a variational 
formulation that is based on extremizing the rest frame kinetic energy of 
flow for fixed relative enstrophy [T2], jlOj, (TH]. This approach is the dual of 
the Minimum Enstrophy method where enstrophy is extremized under fixed 
kinetic energy without fixing angular momentum. The Minimum Enstrophy 
Method is in turn related to the dynamical asymptotic principle of Selective 
Decay for damped 2D turbulence which states that the quotient of enstrophy 
to energy tends to a minimum in time. Amongst many others, the Minimum 
Enstrophy methods with fixed angular momentum were first used by Leith 
28J, Young [2| and Prieto and Schubert in geophysical flows. The main 
result in [T3] gives precise necessary and sufficient conditions for the exis- 
tence and nonlinear stability of prograde and retrograde solid-body flows in 
the coupled barotropic fluid - rotating solid sphere system in terms of the 
planetary spin, relative enstrophy and rest frame kinetic energy. These con- 
ditions have been useful also as starting points in Monte-Carlo simulations 
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of the spherical model for barotropic flows in ^H] and will be compared to 
the results of the mean field theory here as well as the results in [T7j . 

Results from all three methods are closely correlated with recent rigorous 
results on the nonlinear stability of steady-states of barotropic flows [10] , [J2] , 
15J and results from direct numerical simulations of a pde closely related to 
the coupled barotropic fluid - rotating solid sphere system |S]. The rela- 
tionship between nonlinear stability and statisical equilibrium properties of 
geophysical flows is however, not new. Carnevale and Frederiksen discussed 
this relationship in their paper [35], and Shepherd et al |2T] . [22], [2H] are 
good references on applications of the energy- Casimir method in geophysical 
flows. In particular, Shepherd et al [T] , [20] have found convincing arguments 
for zonal anisotropy in rapidly rotating barotropic flows (and the lack of such 
in slowly or non-rotating flows), in the sense that, the Arnold stability theo- 
rems used to prove the nonlinear stability of zonal basic flows, fail when the 
planetary spin is small. 

In conclusion and comparison with the extensive Monte-Carlo simulation 
results on the same problems (cf. Ding and Lim [16J), the numerical asymp- 
totic results for the dissipative coupled barotropic fluid - rotating solid sphere 
system in Cho and Polvani [5 and the variational results in Lim ^3] , we ob- 
serve that (1) the mixed state v = m is preferred for all positive temperature 
except when the planetary spin is large compared to relative enstrophy, whence 
the counter-rotating unmixed state v = nd arises at low enough positive tem- 
peratures, (2) the mixed state v = m is preferred for negative temperatures 
that are more negative than T~ < (when the system has intermediate ki- 
netic energies) while the pro-rotating state v = n u is preferred for negative 
temperatures hotter or less negative than T~ (when the system has extremely 
high kinetic energies), (3a) the positive critical temperature T c increases with 
planetary spin Q and (3b) the negative critical temperature T~ decreases lin- 
early in the planetary spin Q. Its dependence on relative enstrophy is more 
complex, being a sum of two terms, one of which is linear in Q r and other is 
proportional to y/Q^; an increase in Q r results in a decrease in T~ , that is, 
T~ becomes more negative and less hot. 

Given the accepted temperature convention that the more negative T < 
is less hot (or less energetic), and negative temperatures are hotter (or more 
energetic) than positive ones, we note that the Monte-Carlo results of Ding 
and Lim ^H] > the dynamic simulation results of Cho and Polvani [3] and the 
variational steady state results of Lim ^3] are all consistent with the mean 
field results in this paper. Moreover, they complement each other. 
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In particular, when T is negative and hotter than T~(Q, Q r ) (or have rel- 
atively high kinetic energy), the prograde solid-body state v = n u is preferred 
over both the mixed state v = m and the retrograde solid-body state v = n<j, 
according to result (2) in our mean field theory; according to simulations JH] 
of the statistical equilibrium at very hot negative temperatures, the prograde 
solid-body flow is the most probable vortex state; according to the numerical 
results on the dissipative barotropic flows in Cho (cf. fig. 8 in [Hj), there 
is a robust relaxation to the prograde solid-body state when the planetary 
spin Q is not too large for a given range of flow kinetic energy, which is con- 
sistent with result (3b) that T~ decreases linearly in the planetary spin Q; 
and finally according to the steady-state / nonlinear stability results in Lim 
[To] , the prograde solid-body state is allowed only when the planetary spin 
is small compared to the kinetic energy, and in this case, it is nonlinearly 
stable (and hence observable). 

According to result (1) in our mean field theory, the mixed vortex state 
v — m is preferred for all positive temperatures, unless planetary spin is 
large compared to relative enstrophy, in which case, the retrograde solid- 
body state v = n d is preferred over both the mixed state and the prograde 
solid-body state v = n u when positive T is lower than T c (Q,Q r ). According 
to the MC-simulator results in Ding and Lim > the statistical equilibria of 
the barotropic flows relaxes to the retrograde solid-body state for relatively 
low positive temperatures and large planetary spins; this is consistent with 
our result (3a) that T c increases with the planetary spin Q. According to the 
numerical results on the dissipative coupled barotropic fluid - rotating solid 
sphere system in Cho (cf. fig. 9 in j^J), the asymptotic state reached when 
planetary spin is relatively large, is dominated by large anti-cyclonic polar 
vortices; allowing for the fact that the large polar vortex can be expressed 
in the form of a superposition of mainly the retrograde solid-body state and 
some zonally symmetric spherical harmonics of low order, their numerical 
result is partially consistent with our results (1) and (3a). Finally, according 
to the steady-state / nonlinear stability results in Lim jT3], the retrograde 
solid-body state arises for all values of the planetary spin and relative en- 
strophy, but it is nonlinearly stable (and hence observable) only when the 
planetary spin is large compared to relative enstrophy; this is consistent with 
our mean field results. 
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8 Appendix: Non-divergent vs divergent flows 



The non-divergent case of the coupled barotropic fluid - rotating solid sphere 
system in this paper is much simpler to treat than the divergent case which 
is of course more realistic. This more realistic model is essentially based on 
coupling the shallow-water equations (coupledSWE) model to a rotating solid 
sphere through a complex torque mechanism that allows the divergent fluid 
in this case to exchange energy as well as angular momentum with the sphere. 
We will summarize aspects of the coupledSWE in order to state explicitly 
the approximations in the coupled barotropic fluid - rotating solid sphere 
system. For this purpose let us denote by U, L and H, the velocity, length 
and depth scales respectively. Then two important dimensionless numbers 
are the Rossby and Froude numbers respectively, 

r u - u 



2QV vY/T 

where g is the gravitational constant. Within the coupled SWE model, the 
relative importance of convective / local flow to rotational effects is measured 
by the Rossby number R. The Froude number F measures the relative im- 
portance within the coupled SWE model of the convective / local flow effects 
to gravity-depth effects. In a definite sense, a small Rossby number R <C 1 
signals the importance of rotational effects: Q has to be relatively large or 
the scale L of the flow has to be relatively large in order for rotation of the 
planet to be important. On the other hand, a large Froude number F»l 
implies the importance of gravity effects, since in this case, the gravity waves 
are relatively slow, and cannot be time-averaged out of the problem. 
The Rossby Radius of deformation, 

~U RL 



LR ~^m f ' 

measures the relative importance of gravity-depth effects to rotational effects. 
When it is of 0(1), both gravity and rotational effects are relevant to the 
problem, and only when Lr ^> 1 that rotational and convective or inertial 
effects dominate. It is convenient to label the square of Lr/ L by the Burgers 
number 

F 2 L 2 ' 

Small values of B signals the importance of gravity-depth effects over rota- 
tional effects; it inludes the case when the Rossby number R itself is relatively 
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small, that is, when rotational effects dominate convective or inertial effects, 
as well as the case when R is 0(1), that is, when rotational effects are com- 
parable to convective or inertial effects. 

The coupled barotropic fluid - rotating solid sphere system can be charac- 
terized as the limit of the coupled SWE when L R tends to oo or equivalently 
when the depth scale H tends to oo with Q and g fixed. The flow (u, v) in 
the coupled barotropic fluid - rotating solid sphere system model is strictly 
2D, that is, uj is a scalar field and the top and bottom boundary conditions 
are idealized away by taking, in effect, the depth scale if to oo. Thus, in 
a definite sense, the coupled barotropic fluid - rotating solid sphere system 
models a rotating fluid of infinite depth. This fact partly accounts for its 
tractability relative to the more complex coupled SWE model where bound- 
ary conditions at the top and bottom of the fluid are retained. In this definite 
sense, the coupled barotropic fluid - rotating solid sphere system model is 
non- divergent, i.e., div (u,v) = 0, while the coupled SWE is a divergent 
model. 
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